Anomalously localized states and multifractal correlation of critical wavefunctions 
in two-dimensional electron systems with spin-orbital interactions 
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Anomalously localized states (ALS) at the critical point of the Anderson transition are studied for 
the SU(2) model belonging to the two-dimensional symplectic class. Giving a quantitative definition 
of ALS to clarify statistical properties of them, the system-size dependence of a probability to find 
ALS at criticality is presented. It is found that the probability increases with the system size 
and ALS exist with a finite probability even in an infinite critical system, though the typical critical 
states are kept to be multifractal. This fact implies that ALS should be eliminated from an ensemble 
of critical states when studying critical properties from distributions of critical quantities. As a 
demonstration of the effect of ALS to critical properties, we show that the distribution function of the 
correlation dimension D2 of critical wavefunctions becomes a delta function in the thermodynamic 
limit only if ALS are eliminated. 

PACS numbers: 71.30.+h, 72.15.Rn, 64.60.Ak, 71.70.Ej 



I. INTRODUCTION 

The disorder induced metal-insulator transition, 
namely, the Anderson transition^ corresponds to a fixed 
point in a real-space renormalization transformation. 
Approaching the critical point from an insulating phase, 
the localization length £ increases and diverges at the 
critical point. Since the length £ is a unique characteristic 
length near the transition point, there is no length scale 
characterizing the critical state. This implies that the 
wavefunction at the Anderson transition point (the crit- 
ical wavefunction) is scale invariant. In fact, Aokfi and 
Wegner— have demonstrated that critical wavefunctions 
have multifractal properties which can be described by an 
infinite set of exponents. These exponents are quite im- 
portant because they define not only spatial distribution 
of wavefunction amplitudes but also dynamical proper- 
ties of the electron system at criticality^ 

Among various works in which exponents characteriz- 
ing multifractality of critical wavefunctions have been ex- 
tensively s t u died)^iSiL2i2ii24iiii2ii^ numerical multifrac- 
tal analyses based on the box-counting method give the 
most reliable values of exponents at the present stage, 
except for a few exactly-solvable problems^ However, 
values of exponents reported so far widely fluctuate. For 
example, the correlation dimension D% for the three- 
dimensional Anderson model takes values ranging from 
1.45i£to 1.68 A These fluctuations result mainly from the 
following reasons. One is a finite-size effect which can be 
excluded in principle by a finite-size scaling. Another rea- 
son is related to anomalously localized states (ALS)^i& 
at the critical point. Statistical properties of ALS in 
a metallic phase have been studied welli 17 ' 18 i 19 ' 20 i 21 i 22 i 23 
Specific local configurations of disorder in a finite spa- 
tial region may lead to electron states localized within 
this region, while the average disorder is weak compared 
to the critical disorder. Although such ALS induced by 
statistical fluctuations of disorder configurations remain 
even if the system is infinite, they do not contribute to 



ensemble-averaged properties of amplitude distributions 
of wavefunctions or transport properties in the metallic 
phase. This is because ALS have point-like spectra, while 
metallic (extended) states have continuous ones. 

ALS also appear at the critical pointiiiSiS The lo- 
calization length of a multifractal wavefunction diverges 
in an infinite system at the critical point, while sizes of 
ALS are finite. In contrast to the case of metallic phases, 
it is not obvious that ALS can be neglected in infinite 
systems at the critical point. The effect of ALS is much 
stronger than that in metallic phases. However, funda- 
mental knowledge on, e.g., a probability to find ALS or 
influences of ALS to critical properties, are still unclear 
because of the lack of a quantitative definition of ALS. 
It is important to clarify how strong ALS affect critical 
properties by defining ALS quantitatively and revealing 
their fundamental features. 

In the present paper, we give a quantitative and prac- 
tical definition of ALS at criticality suitable for numer- 
ical investigations. In order to define ALS, it is natural 
to evaluate spatial extend of wavefunctions. Quantum 
states with localization lengths less than the system size 
are regarded as ALS in a finite system. This definition 
is, however, not appropriate for precise and systematic 
calculations, because we cannot determine which wave- 
functions are more multifractal if localization lengths of 
two wavefunctions are both close to the system size. We 
show that the correlation function of box-measures of 
multifractal critical wavefunctions&SSiS! works quite well 
for the definition of ALS. In order to demonstrate how 
ALS influence critical properties of the Anderson tran- 
sition, we examine the distribution function of the cor- 
relation dimension D2 of critical wavefunctions. Results 
presented here urge us to reconsider our picture of ALS. 

The efficiency of our definition of ALS is demon- 
strated by applying them to critical wavefunctions in two- 
dimensional electron systems with spin-orbit interactions 
(symplectic systems), which exhibit the Anderson metal- 
insulator transition. For ensuring the multifractality of 
typical critical wavefunctions even in a short-range scale, 
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we adopt the SU(2) model^ in which scaling corrections 
due to irrelevant scaling variables are negligible. 

This paper is organized as follows. In Sec. [HJ we pro- 
pose a quantitative definition of ALS based on the cor- 
relation function of box-measures of wavefunction am- 
plitudes. Some quantities and exponents appearing in 
the multifractal analysis are defined in this section. In 
Sec. IIIII we briefly explain the SU(2) model for which 
the efficiency of our methods is demonstrated and glance 
a numerical technique to obtain eigenstates of the SU(2) 
model. Several ALS calculated numerically for the SU(2) 
model are shown in Sec. II VI We display how well our def- 
inition of ALS works via concrete examples of ALS in the 
SU(2) model. In Sec. El we investigate the fluctuation of 
values of the correlation dimension D% in infinite SU(2) 
systems as an example of the effect of ALS. Section IVll 
is devoted to conclusions. 



II. DEFINITION OF ANOMALOUSLY 
LOCALIZED STATES 

In this section, we propose a quantitative definition of 
ALS based on a multifractal analysis. For this purpose, 
we give some definitions of basic quantities and exponents 
used in the multifractal analysis. At first, we introduce 
a quantity Z(q) defined by 
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where n m = Eieb(i) l^| 2 [or Vb(i) = J m \ip(r)\ 2 dr in 
continuous systems] is the box measure of wavefunction 
amplitudes. The summation Yliebfl) ( or the hitegral 
J b ^ dr) is taken over sites i (or the spatial region) within 
a small box b of size I, and the summation in Eq. (JIJ is 
taken over such boxes. The quantity q is the order of the 
moment Z q (l). For a multifractal wavefunction, Z q {l) 
obeys a power law 



(2) 



where r(q) is called the mass exponent. From Eq. J2J, 
the mass exponent is given by 



t(q) = lim 
/->o 



InZqjV) 
InZ ' 



(3) 



The correlation function G q (l,L,r) defined below can 
describe the box-measure correlations: 8 ' 29 



G q (l,L,r) 



1 



N b N b , 



(4) 



where iJLb T (i) is the box measure of a box b r (l) of size I fixed 
distance r — I away from the box b(l), A& (or Nb r ) is the 
number of boxes b(l) [or b r (l)}, and the summation ^ b is 
taken over all such boxes b r (l). The correlation function 
G q (l,L,r) is a generalized quantity of the qth moment 
Z q (l). For a multifractal wavefunction, G q (l, L, r) should 
behave as£ 



G q (l,L,r) oc l<i)L-y^ r - z{q \ 



(5) 



where x(q), y(q), and z(q) are exponents in multifrac- 
tal correlations. This proportionality relation is sensitive 
to ALS and then suitable for defining ALS. Since our 
purpose is to examine multifractality of individual wave- 
functions by using Eq. JSJ), the system size L is always 
fixed. To find the I and r dependences of G q (l, L,r), we 
concentrate on the following functions, 



J (l) = G 2 (hL,r = l)<xl x W- z W, 



and 



R L {r) = G 2 (l = l,L,r) oc r 



-2(2) 
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In order to quantify non-multifractality of a specific 
wavefunctions, it is convenient to introduce variances 
Var(log 10 QL) and Var(log 10 Rl) from the linear func- 



tions of log 10 1 and log 



ilO ' 



logi 



no QUI) = [x{2) - 
z ( 2 )] lo gio l + c Q and lo gio R U r ) = - z ( 2 ) \og 10 r + c R , re- 
spectively, calculated by the least-square fit. From these 
variances, a quantity T is defined by 

r(L, A) = AVar(log 10 Q L ) + Var(log 10 R L ), (8) 

where A is a factor to compensate the difference between 
average values of Var(log 10 Ql) and Var(log 10 Rl)- Us- 
ing r given by Eq. © , the quantitative and expcdicntial 
definition of ALS at criticality is presented by 



r > r*, 



(9) 



where T* is a critcrial value of T to distinguish ALS from 
multifractal states and chosen appropriately as demon- 
strated later. 



Since ALS are not multifractal, it seems possible to 
define ALS by calculating deviations of the moment Z q (l) 
from the power law Eq. J5J). However, the moment Z q (l) 
represents simply the average of |jL [except for a trifling 

prefactor (L/l) d ], and does not describe the multifractal 
correlation of box measures. Thus, the quantity Z q (l) is 
not sensitive to ALS and Eq. is insufficient for the 
definition of ALS. This will be demonstrated by showing 
an example in Sec. IIVI 



III. MODEL 

In this paper, the efficiency of the above definition of 
ALS and the effect of ALS will be demonstrated for an en- 
semble of wavefunctions at the Anderson metal-insulator 
transition. For this purpose, the effect of scaling cor- 
rections due to irrelevant fields which is a finite-size ef- 
fect and independent of the ALS should be reduced as 
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much as possible. Therefore, considering the advantage 
of system sizes, we focus our attention on the Ander- 
son transition in two-dimensional electron systems with 
strong spin-orbit interactions, in which systems have no 
spin-rotational symmetry but have the time-reversal one. 
Hamiltonians describing these systems belong to the sym- 
plectic ensemble. Several models2fli2ii2S have been pro- 
posed so far to represent the symplectic systems. Al- 
though the Ando model^i has been most extensively 
studied, relatively large scaling corrections make it dif- 
ficult to distinguish the ALS effect. Recently, Asada et 
alM proposed the SU(2) model belonging to the sym- 
plectic class, for which scaling corrections are negligibly 
small. While the SU(2) model is rather mathematical, 
we believe that presented results are qualitatively the 
same with those for other systems exhibiting the Ander- 
son transition. 

The Hamiltonian of the SU(2) model is compactly writ- 
ten in a quaternion representation as 

h = J2 £ 4 c i > ( 10 ) 

i i,j 

where cj (cj) is the creation (annihilation) operator act- 
ing on a quaternion state vector, R^ is the quaternion- 
real hopping matrix element between the sites i and j, 
and Si denotes the on-site random potential distributed 
uniformly in the interval [—W/2,W/2], The strength of 
the hopping V is taken to be the unit of energy. Here- 
after, we denote quaternion-real quantities by bold sym- 
bols. A quaternion-real number x can be written in the 
form 

3 

x = J2 x ^, (11) 

where X IS cL real number and the primitive elements 
r M (^ = 0, 1, 2, 3) define the quaternion algebra as 

t°t° = - T n T n = T° (n = l,2,3), (12) 
rV l = T n r° = T n (n = l,2,3), (13) 

T l T m = -T m T l = T n (14) 

(I, m, n = 1,2,3 and any cyclic permutation). 

In the Hamiltonian Eq. (|10f) . the matrix element Rij is 
defined by 

Rij — cos ctij cos Pij t° + sin 7^ sin /3yT 

— cos ^ij sin fajT 2 + sin otij cos flijT 3 , (15) 

where oiy and 7^ are distributed uniformly in the range 
of [0, 27r), and /3,j is distributed according to the proba- 
bility density P{(3)d(3 = sin(2/3)d/3 for < (5 < it/2. 

We have calculated eigenvectors of the Hamiltonian 
Eq. ljTU|) by using the forced oscillator method (FOM)iS 3 . 
The FOM is an efficient algorism to solve quickly an 
eigenvalue problem with saving memory spaces, in which 



an eigenvector belonging to the eigenenergy E is ex- 
tracted as a resonant mode of the classical dynamical 
system described by the same Hamiltonian to an external 
periodic force with the angular frequency \[E. The FOM 
can be extended to the eigenvalue problem of quaternion- 
real matrices. The corresponding dynamical system is 
described by 

d 2 

— Xi(t) = - HijXjit) + Fi cos(fif), (16) 

3 

where JTy is the quaternion-real matrix element of 
the Hamiltonian Eq. Hl()|) . Xi(t) is the displacement 
vector, and F^ is the site-dependent amplitude of a 
quaternion-real external force. Decomposing the dis- 
placement Xi(t) into quaternion normal modes e^(A) 
as Xi(t) = J2\ e iWQ\(t) 7 the quaternion coefficient 
Q\{t) with A for which E\ is the closest eigenenergy to 
Q? enlarges compared to other coefficients after a long 
time. We should note the order of normal modes e^(A) 
and Q\(t) in the decomposition of Xi(t) reflecting the 
non-commutative property of quaternion numbers. An 
iterative procedure can accelerate the enhancement of 
Q x (t)ml Finally, we can obtain the eigenvector of the 
quaternion-real matrix H as £Ej(i) ~ 6j(A). The obtained 
quaternion-real eigenvector represents two physical states 
simultaneously, which correspond to the Kramers dou- 
blet belonging to E\. Denoting e^(A) by X^ e f( / M rM ' 
the ith elements of the two complex vectors |/(A)) and 
|<?(A)) corresponding to the Kramers doublet are given 
by hi-i{\) = e°(A) + ieJ(A), f 2i (X) = -e|(A) + ief(A), 
Hx-iW = e?(A) + ie\ (A), and g 2i (X) = e°(A) - ie}(\), 
respectively. 34 As we expected, \g(X)) is the time-reversal 
vector of |/(A)), namely, |<?(A)) = — ia y \f(X)), where a y 
is the y component of the Pauli matrix. 



IV. ANOMALOUSLY LOCALIZED STATES FOR 
THE SU(2) MODEL AT THE CRITICAL POINT 

In this section, we demonstrate by showing critical 
wavefunctions of the SU(2) model that ALS can be dis- 
criminated by examining the box-measure correlation 
function defined by Eq. (0J. Critical wavefunctions are 
calculated by the FOM for the SU(2) model with W = 
5.952 for which the eigenstate with E = 1 is known to 
be critical^ All wavefunctions calculated in this paper 
have their eigenenergies closest to E = 1. This implies 
that we extract only one critical state from a system. 
In numerical calculations, the periodic boundary condi- 
tions are employed for both directions. Figure ^ shows 
squared amplitudes of critical wavefunctions in systems 
of L = 120 with different realizations of on-site random 
potentials. The wavefunctions shown in Fig. ^a) seems 
to be strongly localized with a single peak of large am- 
plitudes even for the critical conditions (W = 5.952 and 
Em 1). We can regard this eigenstate as an ALS. In 
order to exclude a possibility that the critical energy is 
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(a) 



(b) 




FIG. 1: Squared amplitudes of critical wavefunctions for the SU(2) model of the system size L — 120. The eigenenergies 
of these states are (a) E — 0.99976, (b) 1.00020, and (c) 0.99970. Wavefunctions labelled by (a) and (b) are spatially localized 
and regarded as ALS. The wavefunction (c) seems to be multifractal. 
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FIG. 2: The functions Ql(1) (filled circles with the left vertical axis) and Rl(t) (open circles with the right axis) for the 
wavefunctions shown in Fig. 1. The labels (a), (b), and (c) correspond to Figs. 1(a), 1(b), and 1(c), respectively. Dashed lines 
are just guide to the eye. 



shifted due to finite-size effects, it has been confirmed 
that eigenstates belonging to adjacent energy levels are 
not localized. The wavefunction shown in Fig. ^b) is 
also localized with a relatively long localization length. 
In this case, we see two peaks of large amplitudes in the 
spatial distribution of squared amplitudes. On the con- 
trary, the wavefunction shown in Fig. ^c) seems to be 
multifractal. The existence of localized wavefunctions at 
criticality are not peculiar to the SU(2) model. We have 
obtained localized wavefunctions also for the Ando model 
at criticality. 

In order to distinguish systematically ALS from wave- 
functions shown in Fig. ^ by the box-measure correla- 
tion function, we calculate QlQ>) and Rl{t) for these 
three wavefunctions. Results are shown in Fig.Oby filled 
[Ql(1)] and open circles [iJi(r)]. Figures Ufa), Hfb), and 
I2tc) correspond to the wavefunctions in Figs. QJa), ^b), 
and ^c), respectively. Dashed lines are guide to the 
eye. The quantity Ql(1) gives the averaged value over 



different ways to divide the system into small boxes of 
size It is found that Ql(1) and Rl{t) for the wave- 
function Fig. ^a) do not obey power laws as shown in 
Fig.[2a), which implies that the wavefunction is not mul- 
tifractal. The function i?i(r) shown in Fig. \2£b) does 
not follow a power law, while Ql(1) seems to be fit to 
a straight line. This is because the wavefunction shown 
in Fig.QJb) has two peaks of large amplitudes and bears 
resemblance to the multifractal wavefunction [Fig. ^c)] 
in a local view. There exists, however, a distinct char- 
acteristic length scale, namely, a distance between two 
peaks. This fact prevents i?i(r) from obeying a power 
law. Since Ql(1) is the same with N^Z 4 (l), Fig.EIb) 
shows that Z q (l) is not sensitive to ALS. On the contrary, 
both Ql(1) and Rl{t) depicted in Fig. |2Jc) well follow 
power laws. This means that the wavefunction shown in 
Fig. ^c) is multifractal as we expected. 

We calculate the distribution function of the quantity 
r defined by Eq. JHJ for ensembles of critical wavefunc- 
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FIG. 3: Distribution functions F(F, L) versus Y for system 
sizes L = 12, 18, 24, 30, 36, 48, 60, 72, 96, and 120, where 
r defined by Eq. JSJ with A = 3 represents the degree of 
non-multifractality of a wavefunction. The inset shows that 
the distribution function F(F, L) is characterized only by the 
average of F. 



tions. We prepare ten ensembles for different system sizes 
(L = 12 to 120). Each ensemble contains 10 4 critical 
wavefunctions. In the definition of T [Eq. ©], we choose 
A = 3, because the average of Var(log 10 Rl) is about 
three times larger than that of Var(log 10 Ql)- Figure 
13 shows the distribution functions -F(r, L) versus V for 
various system sizes L. Recall that eigenstates having 
large values of T are regarded as ALS. As depicted in 
Fig. 3, the function -F(r, L) has a peak at T — for any 
L. This gives an evidence that typical critical states are 
multifractal. A remarkable feature is that the function 
F(T, L) becomes broad as the system size increases. This 
implies that the probability to find ALS increases with L, 
which obliges us to reconsider the role of ALS because the 
influence of ALS to critical properties in infinite systems 
is stronger than that in finite systems. These results 
lead to crucial conclusions as follows: When we study 
critical properties via the average value of quantities X 
defined near the critical point, the average must be the 
geometric mean or the mode value of X, otherwise ALS 
disturb correct information on critical properties. This 
condition has been satisfied by most of previous works. 
Furthermore, when studying critical properties via distri- 
butions of critical quantities such as the level statistics, 
ALS should be eliminated from an ensemble of critical 
wavefunctions to obtain precise information on critical- 
ity. 

We also found that the distribution function F(T, L) 
is characterized only by the (arithmetic) average value 
of r as shown in the inset of Fig. 3. Thus, the function 
F(T, L) can be written in the form: 



F(T,L)=F (T)f(T/f) 



(17) 



where T is the average value of T and is a function of 
the system size L, f is a function independent of L, and 



120.0 
100.0 
80.0 
I (_, 60.0 
40.0 
20.0 



0.0 I 

0.00 



0.02 



0.04 

1/L 



0.06 



0.08 



FIG. 4: System-size dependence of F. Solid line indicates a fit 
by F /Foe = 1 — cL~ v , where Too is the value of F at L — > oo. 



Fq is a normalization factor. The L dependence of f is 
shown in Fig. 4. We see from Fig. 4 that f converges 
to a finite value Too for L — > oo, which implies that the 
probability to find ALS in an infinite system remains to 
be finite. The fact that F(T, L) can be scaled only by f 
would be an important statistical property of ALS. These 
results were obtained for the SU(2) model. However, we 
believe that similar qualitative features will be obtained 
for other systems at the Anderson transition point. 



V. FLUCTUATIONS OF THE CORRELATION 
DIMENSION 

In order to demonstrate how ALS affect critical prop- 
erties in infinite systems, we calculate the distribution of 
the correlation dimension Di of critical wavefunctions of 
the SU(2) model. The correlation dimension D2 char- 
acterizes multifractality of critical wavefunctions and is 
known to be related to some exponents describing dy- 
namical properties at the critical points It had been 
widely accepted that the exponent D2 is universal and 
does not depend on specific samples of critical wavefunc- 
tions. The numerical work of Ref. 35, however, claimed 
that the distribution function of D2 has a finite width 
even in the thermodynamic limit. In response to this sur- 
prising result, recently, several author siL^S^i have veri- 
fied universality of the exponent D2 by precise and large- 
scale numerical calculations. Their results exhibit that 
L>2 takes a definite value independent of samples for infi- 
nite systems, which contradicts the result by Ref. 35. In 
the following, we show that the non-universal property 
of D2 is a consequence of the existence of ALS. 

The correlation dimension D2 is given by 



Z 3 (l) cx I 1 



(18) 



where Z 2 (l) is defined by Eq. (JTJ, thus D 2 = r(2). Al- 
though Eq. I|18|) does not hold for ALS, we calculate 
L>2 even for ALS by force by the least-square fit. The 
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FIG. 5: Distribution functions of D2 for the same ensembles 
as for Fig. 3. Widths of the distribution functions become 
narrow with increasing system size from L — 12 (dashed line) 
to L = 120 (solid line). 
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FIG. 6: Standard deviations of the distribution functions 
P(D<z) as a function of L _1 . Open and filled circles indicate 
results for the original ensembles and the refined ensembles, 
respectively. Solid lines show the scaling fits by Eq. (1191 . 



distribution functions of D2 for the same ensembles as 
for Fig. 3 are shown in Fig. 5. The correlation dimen- 
sion fluctuates over samples in finite systems whereas the 
width of the distribution function P(Dz) becomes nar- 
row as the system size increases. This fluctuation results 
from the existence of ALS and some scaling corrections 
due to a finite system-size. The fluctuation of D2 due 
to ALS remains even in infinite systems. In order to 
obtain the distribution function of D2 for typical criti- 
cal wavefunctions, we eliminate ALS defined by T > T* 
from the original ensembles. To this end, we choose the 
value of r* appeared in the definition of ALS [Eq. 
as T* = 0.03. We study the system-size dependence of 
the standard deviation <r(L) of the distribution function 

Figure 6 shows the standard deviations cr{L) for the 
original ensembles and the refined ensembles (i.e., the 
ensembles in which ALS are eliminated). Error bars of 
a estimated by the bootstrap method are smaller than 
their marks. According to a standard scaling analysis 
that the L-dependence of a statistical quantity at the 
critical point is scaled by an irrelevant length, we can 
write cr(L) as 

a(L) = a x + cL-y, (19) 

where y is an irrelevant exponent and (Too is the standard 
deviation for L — > 00. Fitting data for the original en- 
sembles to Eq. lfH?|) . we obtain = 0.032 ± 0.018 and 
y = 0.39 ± 0.12. The value of (Too is positive finite even 
though its error bar is taken into account. More precisely, 
(Too is not zero with the confidence coefficient 93%. This 
implies that D2 fluctuates even in the thermodynamic 
limit as shown by Ref. 35. On the contrary, fitting data 
for the refined ensembles gives CToo = 0.005 ± 0.013 and 
y = 0.42 ± 0.08. In this case, Coo is very close to zero, 
and the error bar of cr^ is larger than its mean value. 
Therefore, it is natural to deduce that (Too is zero. We 



can then conclude that the correlation dimension D2 for 
typical critical wavefunctions does not fluctuate in the 
thermodynamic limit. It should be noted that ALS were 
also eliminated in the numerical work of Ref. 11 which 
supports the non- fluctuating D2 for L — > 00. 

Finally, we consider the exponent y in Eq. I|19|) . In 
the SU(2) model, the scaling correction due to a spin- 
relaxation length is quite small. What is the origin of 
the irrelevant exponent y in this case? We emphasize 
that the value of y is close to d — D 2 within the numeri- 
cal error, where D2 — 1.66±0.05 estimated for the refined 
ensemble for L = 120. The relation y = d — D2 has been 
analytically predicted by Polyakov^ where y is the ex- 
ponent describing the correction related to the quantum 
return probability. Thus, we claim that the exponent 
y describing the system-size dependence of the standard 
deviation a(L) for the refined ensembles is nothing but 
the irrelevant exponent predicted by Polyakov, which is 
independent of microscopic details of systems. 



VI. CONCLUSIONS 

We have studied statistical properties of anomalously 
localized states (ALS) at the Anderson transition point 
by defining ALS quantitatively. In our definition, a wave- 
function is regarded as ALS when the functions Ql{1) 
and Rl{t) which are the special cases of the box-measure 
correlation functions G q (l,L,r) do not follow the power 
laws Eqs. 10 and Q. Applying the definition of ALS to 
ensembles of critical wavefunctions of the SU(2) model 
which describes a two-dimensional electron system with 
strong spin-orbit interactions, it has been revealed that 
the probability to find ALS at criticality increases with 
the system size, while typical states are multifractal. 
This result suggests that ALS should be eliminated from 
an ensemble of critical wavefunctions if we study critical 
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properties from distributions of critical quantities. We 
also found that the distribution function of T is charac- 
terized only by its average value T, where T quantifies 
non-multifractality of a wavefunction. The influence of 
ALS to critical properties in infinite systems has been 
demonstrated by investigating the distribution of the cor- 
relation dimension D 2 ■ While the distribution function of 
D2 has a finite width in the thermodynamic limit if ALS 
are not eliminated from an ensemble of critical wavefunc- 
tions, D 2 takes a definite value for L — > 00 if ALS are 
eliminated. 

The existence of ALS at the critical point gives cru- 
cial influences also for transport properties at the An- 
derson transition point. The value of conductance for 
ALS is small, while it is relatively large for multifractal 
states. The fluctuation of critical conductance is deeply 
related to the ALS distribution such as Fig. 3. Further- 
more, the frequency dependence of ac conductivity is af- 
fected by ALS as well, because strongly localized ALS 
contribute only to high-frequency ac transport. Our nu- 



merical works have been performed for the SU(2) model. 
We believe, however, that properties of ALS are essen- 
tially the same for systems belonging to other universal- 
ity classes which exhibit the Anderson transition. Fur- 
ther quantitative investigations of ALS reveal concrete 
relations between the nature of ALS and physical phe- 
nomena. 
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